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Abstract 



We demonstrate that the polarization fractions of most tree-dominated B VV 
decays can be simply understood by means of kinematics in the heavy-quark or large- 
energy limit. For example, the longitudinal polarization fractions Rl of the — > 
^ I {Dl'^ , D*+ , p+)D*' and B^ [D*+ ,D*^ , p'^)p^ modes increase as the masses of the 

■ mesons D*^ ,0*^^ emitted from the weak vertex decrease. The subleading finite-mass 

or finite-energy corrections modify these simple estimates only slightly. Our predictions for 
the B — > D^^s^D* polarization fractions derived in the perturbative QCD framework, espe- 
cially Rl ~ 1 for B^ j^*Oj^*o governed by nonfactorizable H^-exchange amplitudes, can 
be confronted with future data. For penguin-dominated modes, such as -B — > p{lo)K* , the 
polarization fractions can be understood by the annihilation effect from the {S — P){S+P) 
operators, plus the interference with a small tree amplitude. At last, we comment on the 
various mechanism proposed in the literature to explain the abnormal B (j)K* polar- 
ization data, none of which are satisfactory. 
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1 INTRODUCTION 



The measured polarization fractions in the B (j)K* decays have exhibited an anomaly. Let 
denote the longitudinal, parallel, and perpendicular polarization fractions of a -B VV 
mode, respectively. It is well known from a helicity argument that these fractions for light 
vector mesons V follow the naive counting rules, 

Rl^I- 0{mllm\) , R\\ ~ ~ 0{m\lm\) , (1) 

if the emission topology of diagrams dominates, where (^y) is the B (V") meson mass. That 
is, Eq. (1) holds for tree-dominated modes, such as i? ^ pp. whose longitudinal polarization 
fraction has been observed to be Rj^ ~ 1 [1, 2]. For penguin-dominated modes, Eq. (1) could 
be modified by annihilation contributions to some extent [3, 4]. However, the B — > ^K* 
polarization fractions shown in Table 1 were found to dramatically differ from the naive counting 
rules, and have been considered as a puzzle. 





Mode 


Pol. 


Fraction 


Belle 


Babar 








Rl 


0.49 ±0.13 ±0.05 [5] 


0.46 ±0.12 ±0.03 [2] 








R± 


0.12i°:J^±0.03 [5] 


5° 






Rl 

R± 


0.52 ±0.07 ±0.05 [5] 
0.30 ±0.07 ±0.03 [5] 


0.52 ±0.05 ±0.02 [6] 
0.22 ± 0.05 ± 0.02 [6] 




Mode 


Pol. 


Fraction 


Belle 


Babar 


B+ 






Rl 




0.96^^:^1 ± 0.04 [2] 








Rl 


0.50 ±0.191^7 [7] 


0.79 ±0.08 ±0.04 ±0.02 [8] 



Table 1: Polarization fractions in the penguin-dominated B — > VV decays. 

In this work we shall investigate all the B — > VV polarization data carefully, and understand 
more the above puzzle. We show that the B — > VV modes can be classified into four categories: 
the first category can be easily understood by means of kinematics in the heavy-quark limit. 
Take the B^ {D*^ , D*^ , p'^) D*~ modes as examples. Under the naive factorization assump- 
tion (FA) [9], QCD dynamics in different helicity amplitudes is absorbed into the universal 
Isgur-Wise (IW) function [10]. The polarization fractions are then completely determined by 
the kinematic factors, leading to the observation that Rl increases from 0.5 to 0.9, when the 
vector mesons from the weak vertex change through D*, D* and p. We shall examine sub- 
leading effects by deriving the perturbative QCD (PQCD) [11, 12] factorization formulas for 
the B D'(^g^D* decays up to next-to-leading power in mo* ^/itib, ttld* ^ being the D*^-^ meson 
mass. To this level of accuracy, the various form factors deviate from the IW function, and the 
nonfactorizable contributions appear. It will be demonstrated that the simple kinematic esti- 
mates are robust under these subleading corrections. As a byproduct, we observe that Rl ~ 1 
for B^ —> D*^D*^, governed by nonfactorizable l^-exchange amplitudes, differs from Rl ~ 0.5 
for other B D*D* modes. This PQCD prediction can be confronted with data in the future. 

The second category is understandable via kinematics in the large-energy limit, which con- 
sists of tree-dominated B meson decays into light vector mesons, such as S ^ {p,Lj)p. Their 
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helicity amplitudes can be expressed in terms of various heavy-to-light transition form factors 
under FA, which are related to each other by the large-energy symmetry relations. It turns 
out that only two universal form factors associated with the longitudinally and transversely 
polarized final states are relevant. If these two form factors do not differ much, the polarization 
fractions will be completely determined by the kinematic factors, leading to ~ 1 consis- 
tent with Eq. (1). The same argument applies to the {D*^ , D*^)p'^ decays with the 
mesons being emitted from the weak vertex, whose Rl ~ 0.7 are predicted, and can 
be compared with the future data. We then examine subleading corrections to the large-energy 
symmetry relations by including the two-parton twist-4 contribution. Adding this piece to 
the longitudinal polarization amplitude makes complete the ncxt-to- leading-power analysis at 
two-parton level, since the transverse polarization amplitudes are of next-to-leading power by 
themselves. It will be demonstrated that this subleading effect is also neghgible. 

The third category contains penguin-dominated modes, such as S — > p{uj)K*, which can 
exhibit a sizable deviation from the naive counting rules in Eq. (1). As shown in [3], the 
annihilation amplitudes associated with the {S — P){S -\- P) operators follow the counting rules, 

Rl ~ R\\ ~ . (2) 

The PQCD analysis has indicated that the penguin annihilation contribution, together with 
nonfactorizable corrections, bring the longitudinal polarization fraction in a pure-penguin mode 
from Rl ~ 0.9 down to 0.75 [3]. Therefore, the — * p^K*^ polarization data in Table 1 can 
be well accommodated within the Standard Model, showing no anomaly. The longitudinal 
polarization fraction of another mode — > p^K*^ remains as Rl ~ 0.9 because of the inter- 
ference between the penguin amplitude and an additional tree amplitude. From the viewpoint 
of PQCD, the B — > ujK* decays do not differ much from B — > p^K*, and are expected to show 
similar Rl- This prediction can be tested by the future data. 

The fourth category, involving the B (j)K* decays, is the abnormal one. These decays 
occur mainly through the penguin operators, but their Rl are as small as 0.5, much lower 
than 0.75 expected from PQCD. The mechanism proposed in the literature to explain the 
B — > (f)K* polarization data includes new physics [13], the annihilation contribution [14] in 
the framework based on the QCD- improved factorization (QCDF) [15], the charming penguin 
[16], the rescattering effect [17. 18, 19], and the h — ^ sg transition (the magnetic penguin) [20]. 
We shall comment on these proposals: the annihilation contribution has to be parameterized 
in QCDF, and varying free parameters to fit the data can not be conclusive. The charming 
penguin strategy, demanding many free parameters, does not help understand dynamics. The 
rescattering effect is based on a model-dependent analysis [21, 22], and constrained by the 
B — > pK* data. The magnetic penguin is suppressed by the G-parity, and not sufficient to 
reduce Rl down to 0.5. Therefore, none of these proposals is satisfactory [4]. However, we are 
not concluding that the B — > (pK* polarization data signal new physics, since the complicated 
QCD dynamics in S — > VV modes has not yet been fully explored. 

In Sec. II we study the kinematic effects on the polarizations of tree-dominated decays 
using FA in the heavy-quark or large-energy limit. The next-to- leading-power corrections to 
the B D'(^g^D* decays, and the two-parton twist-4 corrections to the B ^ pp decays are 
calculated in Sec. III. We comment on the proposals for explaining the abnormal B — > (pK* 
data in Sec. IV. Section V is the conclusion. 
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2 NAIVE FACTORIZATION 



In this section we demonstrate that the polarization fractions of tree-dominated B — > VV 
decays can be simply understood by means of kinematics in the heavy-quark or large-energy 
limit. 



2.1 Heavy-quark Limit 

We first investigate the polarizations in the decays {D*^ , D*^ , p^) D*^ in the heavy 

quark limit. These modes are color-allowed with the D*'^ , D*~^ , p'^ mesons emitted from the 
weak vertex, respectively, and FA is supposed to work well. Take the S° — > D*'^D*~ decay as 
an example. The B meson momentum Pi, the D* meson momentum P2, and the D* meson 
momentum P3 are chosen, in the light-cone coordinates, as 

Pi = ^(1,1, Or), 

P2 = -^{r2V :r2ri , Ot) , 

Ps = ^(l-r2r;+,l-r2r?-,0T), (3) 

where the factors 77^ are defined by 77''= = rj ± y/rj'^ — 1 with rj = vi ■ V2 being the velocity 
transfer, vi = Pi/uib and V2 = P2/mD* the B meson and D* meson velocities, respectively, and 
r2 = mD*/TaB the mass ratio. To extract the helicity amplitudes, the following parametrization 
for the longitudinal polarization vectors is useful: 

e2{L) = ^2 - = ^(77+, -77-,Ot) , 

P2 ■ n_ V2 

= v,-^^n^^^(--^,l-r,v-,OT] , (4) 
P3 ■ n+ V2r3 \ l-r2ri J 

with the D* meson velocity ^3 = Ps/uid*, the mass ratio — ruD'/mB, and the null vectors 
n+ = (1,0, Ot) and n_ = (0,1, Ot), which satisfy the normalization e|(L) = eK-^) = —1 and 
the orthogonality e2{L) • P2 — €3{L) ■ P^ — 0. For the transverse polarization vectors, we simply 
choose 

e2(T) = (0,0,lT), e3(T) = (0,0,lT). (5) 
The relevant effective weak Hamiltonian is given by 

Kff = ^ Va,V:^[Ci{p.)Oi{p.) + C2(/x)02(/i)] , (6) 

where F's are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, Ci and C2 the Wilson 
coefficients, and 

Oi = {sb)v-A{cc)v-A , O2 = {cb)v-A{sc)v-A , (7) 
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the four-fermion operators with the definition {qiq2)v-A = ?i7/i(l — 75) ?2- The contributions 
from Oi and O2 can be combined, and the resultant coefficient appears as ai — C2 + Ci/Nc, 
Nc being the number of colors. 

The — > D*'^D*~ decay amphtude in FA is expressed as 

^ (D;+(P3,e*)D*-(P2,e;)|Kff|i?°(Pi)), 

x(L>*-(P2,e;)|67'^(l-75)c|S°(P0), (8) 

where the superscript a denotes a possible final helicity state. The first matrix element defines 
the D* meson decay constant, 

{D:^{Ps, el)\c^,{l - 75)s|0) = mn^jDtel . (9) 

The matrix elements for the B ^ D* transitions are parameterized as 

{D*-{P2,e;)\h^^,c\B'{P,)) = V'^?^^^[U{v){v + l)er-U2{vy2-viv^ 

-U3{vy2-viv^] , (10) 

where the form factors ^Ai, ^a^, ^As, and satisfy the relations in the heavy-quark limit, 

Cv = U = Us = C, a = , (11) 

with ^ being the IW function [10]. 

The B'^ D*'^D*~ decay rate is given by 

r = — ^ VA^('^)t;tl('^) , (12) 
87rm^ V 

where Pc = \P2z\ = I-P32I = iTLBr2Vv'^ ~ 1 is the momentum of either of the vector mesons. The 
amplitude Ai^'^^ is decomposed into 

M^^) = {mlML: mlMNeliT) ■ 6^(7), iMre^^^^^'eLe^^Ps^Pap) , (13) 

where the first term corresponds to the configuration with both the vector mesons being longi- 
tudinally polarized, and the second (third) term to the two configurations with both the vector 
mesons being transversely polarized in the parallel (perpendicular) directions. The helicity 
amplitudes are then defined as, 

Al = -GmlML, 



= GmD*mD*j2[{v2-V3y-l]MT, (14) 
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with the normahzation factor G = y Pc/{8TTm%T), which satisfy the relation, 

+ + = (15) 

It is easy to read the hehcity amphtudes off Eqs. (8)-(10), 

Al oc e;{L) ■ el{L){n + 1)U - 4{L) ■ v,[el{L) ■ v,U. + 4{L) ■ V2Us] : 
A|| oc -V2{ri + , 

oc -rs^2[(v2 ■ - l]^v ■ (16) 

In the heavy-quark hmit, i.e., applying Eq. (11), all the hehcity amplitudes depend on 
a single IW function ^, which absorbs QCD dynamics. Simply inserting ttib — 5.28 GeV, 
mo* — 2.11 GeV, and uid* — 2.01 GeV into the kinematic factors, we obtain 

Rl ~ 0.52 , ~ 0.43 , R± ~ 0.05 . (17) 

Equation (16) is apphcable to the 5° D*^D*~ and p^D*~ decays by substituting 

mD* and rUp — 0.77 GeV for mD*, leading to 

Rl ~ 0.54 , R\\ ~ 0.41 , R^ ~ 0.05 , (18) 
Rl ~ 0.88 , R\\ ~ 0.10 , R^ ~ 0.02 , (19) 

respectively. All the above results are consistent with the observed values listed in Table 2, and 
with the predictions in [29]. The estimated polarization fractions for the decay B~ — > K*~D*^ 
are close to Eq. (19) and consistent with the data. For the B^ — > p^D*^ decay, the internal- 
W emission amplitudes involving the B ^ p form factors are suppressed by the vanishing 
coefficient 02 = Ci + C2/NC. Hence, this mode is similar to B^ p~^D*~, and the result 
in Eq. (19) applies as shown in Table 2. It is easy to find that the longitudinal polarization 
fraction increases as the mass of the vector meson from the weak vertex decreases. 



2.2 Large-energy Limit 

We then show that the polarization fractions in the B — >• {p,uj)p decays can be understood by 
means of kinematics in the large-energy limit. The parametrizations of the momenta in Eq. (3) 
and of the polarization vectors in Eqs. (4) and (5) hold here, with the mass ratio r2 for the 
vector meson from the B meson transition and for the vector meson emitted from the weak 
vertex. The transition form factors associated with a B ^ V transition are defined via the 
matrix elements, 

{ViP2,€*2)\h'Q\BiPi)) 
{V{P2,el)\br^,q\B{P,)) 



^'^i^ ^ 6^-^-6LP2,Pi., (20) 

TUB + my 

2mvAo{q'f^q^ + (m^ + mv)A^{q') (^eT - ^g^) 
-A^iq')^^ (Pr + - ^^i^g'^) , (21) 

TUB + my V Q J 
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Mode 


Pol. 


Pr action 


Data 




B+ p+D*^ 
B- K*-D*° 




Rt 
Rl 
R± 
R± 
Rl 
Rl 
Rl 


0.52 ± 0.05 [23] 
0.57 ±0.08 ±0.02 [24] 
0.19 ±0.08 ±0.01 [24] 
0.063 ± 0.055 ± 0.009 [23, 25] 
0.885 ±0.016 ±0.012 [23] 
0.892 ±0.018 ±0.016 [23] 
0.86 ± 0.06 ± 0.03 [26] 




Mode 


Pol. 


Pr action 


Belle 


Babar 


B+ ^ p+pO 
B^ p+p- 
B+ p+u 




Rl 
Rl 
Rl 


0.95 ±0.11 ±0.02 [1] 


0.97l^:;^-^ ± 0.04 [2] 
0.99 ± OMtom [27] 
0.88^2:11 ± 0.03 [28] 



Table 2: Polarization fractions in the tree-dominated B — > VV decays. 



with the momentum q — Pi — P2. The form factors V, Ai, and A2 satisfy the symmetry 
relations in the large-energy limit, 



rriB + my 



2E 



2E 



(22) 
(23) 



where E is the V meson energy, and the function C\\{E) {^±{E)) is associated with the transition 
into a longitudinally (transversely) polarized V meson. Our definition of differs from those 
in [30, 31], such that we have ^± ~ ^y, when not distinguishing the longitudinal and transverse 
polarizations. 

Consider the 5+ p'^p^ mode as an example, which is dominated by a tree contribution 
(assuming that the electroweak penguin contribution is negligible). The explicit expression of 
the relevant effective weak Hamiltonian will not be shown here. In PA, the decay amplitude is 
written as 



(pVl^effl^-') 



ai(p+(P3,6;)|«7;.(l -75)c^|0)(p°(P2,e;)|67'^(l -75)«|5+(Pi)) 



±a2(p°(P2,6;)|iZ7^(l-75)^i|0)(p+(P3,6;)|67'^(l-75)c?|5+(A)) 



(24) 



where the two terms can be combined, leading to the helicity amplitudes, 

2^2 ■ -P3-P2 ■ ^3 



Al oc (l±r2)e;(L)-e;(L) 



m%{l + r2Ye*2 ■ el 



A, 



A\\ oc 
A 



-^(l + r2)A , 



oc 



2r2r3 
l±r2 



p[{v2-v^Y-l\V . 



(25) 
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The relative phases among the hehcity amphtudes arc (f)\\ = 0^ = tt under FA in our convention. 
In the large-energy limit, i.e., employing Eqs. (22) and (23), Eq. (25) becomes 



Al oc 2r2el{L) ■ el{L)i\\ , 
A\\ oc -2v^r2?7C± , 

A^ oc -2r2r^^j2[{v2-v^Y-l\U , (26) 

implying that QCD dynamics has been absorbed into the two functions and For a rough 
estimate, adopting ^ (/p//J)C±, fp = 200 MeV and /J = 160 MeV [32], the B+ p+p^ 
polarization fractions are given by, 

Rl - 0.95 , Rii - 0.03 , R± - 0.02 , (27) 

consistent with the data in Table 2. The polarization fractions of other tree-dominated B — > VV 

modes, inchiding p'^p~ and p^u, can be explained in a similar way. 

We generalize Eq. (26) to the {D*^,D*^)p^ modes, which arc mainly governed by 

the B p form factors, with the masses rriD*, ttid* being substituted for nip, respectively. 
Their polarization fractions are predicted to be, 

D*V: i?L~0.70, i?||~0.16, R±r^O.U, 

D*+p^ : Rl - 0.72 , R^ - 0.15 , R± - 0.13 , (28) 

which can be compared with the data in the future. 

3 SUBLEADING CORRECTIONS 

Away from the heavy quark limit, the form factors in Eq. (10) deviate from the IW function. 
Beyond FA, nonfactorizablc contributions appear. Both corrections, being subleading [33, 34], 
can be calculated more reliably in the PQCD approach based on kx factorization theorem 
[35, 36] due to the stronger end-point suppression in the former [33] and to the soft cancellation 
between a pair of nonfactorizable diagrams in the latter [37]. In this section we shall examine 
whether the simple estimates made in the previous section are robust under these subleading 
corrections. The factorization theorem for the B D* form factors in the large-recoil 
region of the D* meson can be proved following the procedure in [38] , which are expressed as 
the convolution of hard kernels with the B and D* meson wave functions in both the momentum 
fractions x and the transverse momenta kr of partons. A hard kernel, being infrared-finite, 
is calculable in perturbation theory. The B and D* meson wave functions, collecting the soft 
dynamics in the decays, arc not calculable but universal. After including the parton kx, the end- 
point singularities, which usually break QCDF at subleading level, do not appear, and PQCD 
factorization formulas are well-defined. This formalism has been applied to the semileptonic 
decay B D^*Hu [33], and the nonleptonic decays B D^*h{p) [34, 39] and B 
[40] successfully. 
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V -A 
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V - A 
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O 
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(a) (b) (c) (d) 

Figure 1: Lowest-order diagrams for the — > D*^D*^ decay. 



3.1 



D*+D* 



For the D*~^D*~ mode, we shall neglect the penguin contribution, which is suppressed 

by the Wilson coefficients Ca-g ~ O.OSai. According to the lowest-order external-VT emission 
diagrams in Fig. 1, the decay amplitudes for different final helicity states are expressed as 



71 



cb'' cs 



N 



I e 



af3-yp 



(29) 



where Tl,n,t come from the factorizable diagrams. Figs. 1(a) and 1(b), and M.l,n,t from the 
nonfactorizable diagrams. Figs. 1(c) and 1(d). 

We ffist compute each factorizable amplitude in terms of the "form factors" as an expansion 
in r2. 
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iA2 



Am 

SA2 



(NL) 



(30) 



where the superscript NL denotes the next-to-leading-power corrections. Equation (30) is 
equivalent to the heavy-quark expansion of the heavy- heavy currents in 1/mf, and in l/m^ [41], 
nib i^c) being the b (c) quark mass. We refer the explicit expressions of the PQCD factorization 
formulas for ^ and ^f^^'' to [33], which have absorbed the Wilson coefficient oi in the current 
analysis. In terms of Eq. (30), the factorizable amplitudes are written as 



= 

!Fn = 
Tt = 

A numerical analysis gives 

eJ^)~-o.o2e 



r2rs iv + 1) + ^T) , 



.(NL) 



-0.19 e, i 



(NL) 



-0.05 e 



(31) 



(32) 
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Inserting the deviation from the IW function into Eq. (16), the polarization fractions are mod- 
ified only slightly: 



Rl - 0.54 , Ru - 0.41 , Rj 



0.05 . 

.(NL) 



(33) 



The largest next-to-leading-power correction comes from , which is, however, suppressed 

by the factor r2. 

Next we calculate the subleading corrections from the nonfactorizable diagrams in Figs. 1(c) 
and 1(d). For simplicity, we shall expand all the kinematical factors and the polarization 
amplitudes up to next-to-leading power in r2,3. The factorization formulas are given by 



A4l = IGnCF^/^Ncm^ j dxidx2dx3 j 



M 



N 



M, 



bidbi b^dba (j)B {xi , 6i) (j)D* {X2) (t>Di {x^) 

X (x3-xi + r2X2) -Bfe (4^^ ) /ift^^ (2^1,3^2, 3:^3, &1,&3) 

- (1 -FX2 - xa - xi - X2r2) Eb{tf'') h'^\xi,X2,X3,bi,b3) , 

l&T:CF\l'2.Ncm\ dxidx2dx3 63^63 0b(xi, 61) 0£).(a;2) (0:3) 

X xs r3 Eb{ti^^) hi^\xi, ^2, X3, 61, 63) 

- (1-^x3) r2,Eb{tf'') hf\xi,X2,X3,bi,b3) 
2Mn ■ 



(34) 



(35) 
(36) 



Radiative corrections to the meson wave functions generate the double logarithms In^ kx 
from the overlap of coUinear and soft enhancements, whose Sudakov resummation has been 
studied in [42] . The Sudakov factors from kr resummation for the B meson, the D* meson and 
the D* meson are given, according to [43], by 



exp[-5s(/x)] 

exp[-SD*{l^)] 
exp[-5n*(/i)] 



exp 



exp 



exp 



-s(A;i ,61) - ^ / ^7(«.(/^)) 
3 Ji/bi jJi 

3 Ji/b2 

-s(/c3 ,63) ^7(a.(/^)) 



(37) 



respectively, with the quark anomalous dimension 7 = —as/n. The momenta fcf = XiP{', 
^2 = X2P2 ^ and fcg' = a^aPg", carried by the light valence quarks in the 5, D*, and mesons, 
respectively, define the momentum fractions x. The impact parameters 61, 62, and 63 are 
conjugate to the transverse momenta carried by the light valence quarks in the B, D*, and D* 
mesons, respectively. For the explicit expression of the Sudakov exponent s, refer to [11]. Note 
that the coefficient 5/3 of the quark anomalous dimension in Eq. (37) differs from 2 for a light 
meson. The reason is that the rescaled heavy-quark field adopted in the definition of a heavy- 
meson wave function has a self-energy correction different from that of the full heavy-quark 
field [43]. The evolution factors are then given by 



(38) 
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with the Sudakov exponent S = Sb + Sd* + Sd*- 
The functions h[^\ j = I and 2, are written as 

h^b^ = [e{h-bs)Ko{BmBh)Io{BmBb: 



•3) 



+9{bs - bi)Ko (BniBbs) /q {BruBbi)] 
( KoiBjmBbs) forS|>0 

[ fH^'^ (v^^Bfta) for S| < ' ' ^^^^ 



with the variables 

B'^ = XiX2r2r]'^ 



Bl = xiX2r2'n'^ - X2Xz {r2'n'^ - rf^ , 

Bl = xiX2r277+ -X2{l-xz) (r277"^ - r'^ + (xi + xz) (l - r2?7+) + x^ {rl - ra??") . (40) 
The scales t^J^ are chosen as 



4'^ = max (SmB,<y^mB, 1/61,1/63) ■ (41) 
The B and D* meson wave functions involved in Eqs. (34)-(36) are also referred to [33]: 



(j)B{x,b) — Nbx {1 — x) exp 





/ xniB^ 


2 , ,2 l2' 






' 2 



(42) 



M^) - -^x(l-x)[l + Cz,*(l-2x)] . (43) 

The D* meson wave function 0/). will be assumed to have the same functional form as (po*- 
We do not consider the b dependence of (po* in the large recoil region of the D* meson [33]. 
The Gaussian form of 0b was motivated by the oscillator model in [44]. The shape parameter 
lob — 0.40 GeV comes from [45], and Co* — 1-04 is determined from the IW function ^{r] — 
1.3) = 0.7. The normalization constant Nb is related to the decay constant /b through 



/ 



d.M:^.6 = 0) = ^. (44) 



We shall not distinguish the D* meson wave functions associated with the longitudinal and 
transverse polarizations. There are various models of the B meson wave functions available in 
the literature [43]. It has been confirmed that the model in Eq. (42) and the model derived in 
[46] with a different functional form lead to similar numerical results for the B ^ tt form factor 
[47]. 

Taking into account only the nonfactorizable corrections, the polarization fractions become 

Rl - 0.54 , R\\ - 0.40 , R± - 0.06 . (45) 
Including both the subleading factorizable and nonfactorizable corrections, we have 

Rl ~ 0.56 , R\\ ~ 0.39 , Ri_ ~ 0.06 , (46) 
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Figure 2: VT-exchange diagrams for the — > D*^D* decay. 



which are still close to the values in Eq. (17). Hence, the simple kinematic estimate made in 
the heavy-quark limit is very reliable. 

At last, we compute the — >■ D*^D*^ branching ratio, considering only the leading 
contribution. Employing the CKM matrix elements Vd, = 0.0412, Vcs = 0.996, and Vcd = 
—0.224, the quark masses = 4.8 GeV and rrit = 174.3 GeV, the meson decay constants 
fe = 200 MeV, fo- = 230 MeV, and /^j = 240 MeV, the lifetimes tbo = 1.542 x lO^^^ 
and tb± = 1.674 x 10~^^ sec, and the Fermi constant Gp = 1.16639 x 10^^ GeV^^, we predict 



D*+D*- 



(2.6 



+1.1 
-0.8 



% 



(47) 



which is consistent with the recent measurement (1.85 ± 0.09 ± 0.16)% [48]. The theoretical 
uncertainty in Eq. (47) arises from the allowed range of the shape parameter ujb = (0.40 ±0.04) 
GeV [45] in the B meson wave function. Note that the polarization fractions are insensitive to 
this overall source of uncertainty. 



3.2 B D*D* 

For the B^ —>■ D*^D*~ mode, there exists an additional contribution from the VT-exchange 
topology shown in Fig. 2 compared to B^ — > D*~^D*~ . The factorizable py-exchange diagrams. 
Figs. 2(a) and 2(b), vanish exactly because of the helicity suppression. An explicit evaluation 
shows that the nonfactorizable VT-exchange diagrams. Figs. 2(c) and 2(d), contribute only 
2% of the external- ly emission ones in Fig. 1. Therefore, the VT-exchange topology is even 
less important than the penguin one. Note that the VT-exchange contribution is as important 
as the nonfactorizable contribution for the B — Z)'^*^7r(p) decays [34], to which the helicity 
suppression does not apply. The B^ —>■ D*^D*~ factorization formulas are then similar to 
those of B^ D*^D*^ but with the appropriate replacements of the D* meson mass, decay 
constant, and wave function by the D* meson ones, respectively. Similarly, if neglecting the 
l^-exchange and penguin contributions, the B^ —>■ decay amplitudes will be the same 

as of 5° ^ D*+D*-. 

The numerical results are summarized below. Including both the subleading factorizable 
and nonfactorizable corrections, we obtain the polarization fractions of the 5° D*^D*^ and 
5+ D*+D*^ modes, 

Rl ~ 0.56 , R\\ ~ 0.38 , ~ 0.06 , (48) 
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Figure 3: Lowest-order diagrams for the f)*0£)*o jg^ay, 



which are also close to the simple estimate in Eq. (18). Using the same input parameters for 
the quark masses, the B meson lifetimes,..., we predict the branching ratio. 



B(B^ 



1.2iOJ) xlO- 



(49) 



which is consistent with the updated measurement B{B^ D*^D*^) = (8.1 ±0.8 ±0.1) x 10~^ 
[24]. The branching ratio B{B^ D*~^D*^) can be obtained simply by changing the B meson 
lifetime, whose value is close to that in Eq. (49). 

Since only the ly-exchange topology in Figs. 3(a)-(d) and the penguin annihilation topology 
in Figs. 3(e)-(h) contribute, the B^ f)*Of)*o (jgc^y must have a tiny branching ratio. For 
both topologies, the factorizable amplitudes diminish because of the hehcity suppression (only 
the penguin annihilation from the {S — P){S + P) operators survives, which do not exist in 
this mode). Hence, they are not exhibited in Fig. 3. With the penguin amplitudes being 
down by the Wilson coefficients, we shall calculate only the tree contribution. By measuring 
the B^ £)*0£,*o mode, we learn how the nonfactorizable VT-exchange contribution governs 
the polarization fractions. It will be shown that the polarization fractions in this decay differ 
dramatically from those in the _B° D*^D*^ and 5+ /}*+/)*o ^lecays. Qn the other 
hand, it has been proposed [49] to extract the weak phase 03 from the B — ^ D^*^D'^*\ D'^^D'^*^ 
measurements. 

The amplitudes from Figs. 3(a) and 3(b) are given by 



Ml = 167rCir-y/2iVcm| J dxidx2dx^ J 



bidbi h2dh2 <Pb{xi, h) (po- {X2) (pD* (2:3) 



(xs -xi) Ef{tf^)hf\xi,X2,X3,bi,b2) 
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-X2 Ef{tf) hf\xi,X2, X3,bi,b2) 

Mn,t = O (r^) , 
and those from Figs. 3(c) and 3(d) are 



(50) 
(51) 



M'l = 167rCF^/2Ncm% dx-idx2dx3 62(^^2 0B(a:^i, &i) 0d* (2^2) 0d*(2;3) 



X 



(1 -X2) Ef(t}^^) hj^\xi,X2,X3,bi,b2) 



- (1 - X3 + xi) Ef{tP) hj'^\xi, X2, X3, bi, 62 



M'r 



O {r^ 



(52) 
(53) 



It is observed that the transverse polarization amplitudes vanish in the current next-to-leading- 
power accuracy. The evolution factors are 



Ef{t) = «.(t)^exp[-5(i)|,3=,,] . 



(54) 



The hard functions h^P^^\ j — ^ and 2, are written as 



h 



cm 



X 



e{b^ - b2)H^^^ (F^'^mBbi) Jo (FWmB62) 
+eib2 - bi)H^'^ (FWrnB^s) Jo {F^'^msbi) 
( Ko (F^mBbi) for F^^ > 
H^'U^lFf^lrnnb,) for i^^^ < 



VK 

2 



(55) 



with the variables 



F^ 

?2 



X2r2rj^X3{l -r2r) ) , 



= X2r2r] 



p2 
-^2 

77/2 



^,2 



xi - xs(l - r277~) 
1 + X2r2r]^) xi - 1 2:3(1 - r2r}~) 

{1-X2r2r]^) 1 - X3{1 - r2rj~) , 
(1 - X2r2ri~^) X1-I + ^3(1 - r27]~) 
1 + X2r2r)^ xi - X3{1 - r2r]~) 



(56) 



The variables F"s can be obtained from F's by interchanging X2r2V^ and 1 — X2r2r]~^, and 



2^3(1 — r2r) ) and 1 — ^3(1 — r2r) ). The hard scales t^^^^ are chosen 



as 



(^FWms, ^JlFf^lrriB, I/61, 1/62^ ■ 



(57) 
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The longitudinal polarization fraction dominates: 

Rl^I, ~ ~ few% , (58) 

which differs very much from that in Eq. (48). Therefore, the comparison of the theoretical 
prediction with the future data can test the PQCD approach. We also predict the branching 
ratio, 

5(^0 ^ ^*o^*o^) _ (8.9t\i) X 10-^ , (59) 
which can also be compared with the future data. 



3.3 B 



PP 



In this subsection wc examine whether the simple estimate of the polarization fractions in 
the tree-dominated B meson decays into two light vector mesons is robust under subleading 
corrections. Similar to the previous subsection, the 0{mv/mB) terms should be included into 
the factorizable amplitudes at this level of accuracy. At the same time, the two-parton twist-4 
contribution appears, since the linear end-point singularity involved in collinear factorization 
theorem modifies the power behavior from 0(my/m^) into 0{vi\/ /ms) [45]. The inclusion of 
these two corrections makes complete the next-to-leading-power analysis at the two-parton level. 
The nonfactorizable amplitudes have been known to be small due to the strong cancellation 
between a pair of nonfactorizable diagrams [11, 12]. The annihilation amplitudes for tree- 
dominated modes are also negligible due to helicity suppression [37]. Hence, we shall not 
consider the two-parton twist-4 correction to these two subleading contributions. Below we 
analyze the B ^ pp longitudinal polarization amplitude as an example. 

The two-parton p meson distribution amplitudes up to twist 4 are defined by the following 
expansion [32], 



(p-(P2,e;j|d»X0)/|0) = 



1 



•2-Z 



+mpl(t)l{x) - 



i)2(t>p{x)+mp{fi+ /i_ - 1)0* (x) 



2p2 ■ n_ 



(60) 



ij 



where the new vector p2 contains only the plus (large) component of P2. The distribution 
amplitude 4>p is of twist 2 (leading twist), 0^ and (p^ of twist 3, and 0^ of twist 4. The twist-3 
distribution amplitudes in fact give leading-power contribution due to the similar modification 
from the end-point singularity. The explicit expressions of the above p meson distribution 
amphtudes are referred to [45], and 0^ is given by [50] 



1 - 1.62Cl^\2x - 1) - 0AlC^^"{2x - 1 



2V2N, 

with the Gegenbauer polynomials. 



1/2/ 



Cl^'it) = ^{3t' - 1) , C^'it) = ^{35t' - 30t' + 3) . 



(61) 



(62) 
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The longitudinal factorizable amplitude in the B ^ pp decays is written, up to twist 4, as, 

/■I poo 

Tl = SnCFm^ / dxida;2 / bidbib2db2 (pB{xi,bi) 
Jo Jo 

X I [((1 + X2)(l - rl) - (1 + 2x2)rl) 0p(x2) 

+r2(l - 2x2) (0p(a;2) + ^^(a^s))] £;e(^i'^)/ie(^i, ^2, 61, 62) 
20^(0:2) +r2 0^(a;2)l E,{tf^)h,{x2, x,,b2,b^)] , 



+r2 



(63) 



where the first (second) term containing Eg (4^'') [-^e(4^'*)] comes from the lowest-order diagram 
similar to Fig. 1(a) [Fig. 1(b)], but with the D* and D* mesons being replaced by the p mesons. 
The evolution factor is 



E,{t) = a,{t) ai(t) exp[-SB{t) - Sp{t)] , 
with the Sudakov factor from the kr resummation, 

dp, 



ex-p[-Sp{iJ,)] ^ exp 



-s{k+, 62) - s{P,^ - kt, 62) - 2 r 4l{aM) 

Jl/b2 jJ. 



and the hard scales, 

= m.&x.{^/x^mB, 1/61, 1/62) , tf^ = m.six.{^/xlmB, I/61, 1/62) 
The hard function is given by 



(64) 



(65) 



(66) 



he{xi,X2, bi, 62) = St{x2) Kq {^/x^niBbi) [9{bi - b2)Ko {^/x^rriBbi) Iq {^/x^mBb2) 

+9{b2 - bi)Ko {y/x^mBb2) Iq (y/x^mBbi)] . (67) 

The Sudakov factor St{x) arises from the threshold resummation of the double logarithms 
as In^ X, which are produced by the radiative corrections to the hard kernels. Its expression 
[51], 



St{x) 



2i+^T(3/2 + c) 
0Fr(l + c) 



[x{l-x)]\ 



(68) 



with the constant c ~ 0.3, provides further suppression in the end-point region of a; — 0, 
and improves the perturbative calculation. Since we have performed the leading-logarithm 
resummation so far, only the behavior of St{x) at small x is reliable. The parametrization in 
the large x region, also vanishing, was proposed for convenience [51]. 

The numerical results are listed in Table 3, where diagram (a) [diagram (b)] refers to the 
lowest-order diagram with the hard gluon being on the B (p) meson side. It indicates that 
the next-to-leading-power corrections decrease (increase) the leading-power contribution from 
diagram (a) [diagram (b)] slightly. Considering the net effect, these corrections are indeed 
neghgible. 
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Diagram (a) 


Diagram (b) 


sum 


Leading-power 


0.330 


0.088 


0.418 


Plus next-to-leading power 


0.324 


0.096 


0.420 



Table 3: Contributions to Tl up to next-to-leading power. 



4 COMMENTS ON PLAUSIBLE EXPLANATIONS 

We now comment on the mechanism proposed in the literature to explain the abnormal polar- 
ization fractions of the B — > (^K* decays. It will be argued that these proposals involve many 
free parameters, can not account for the polarizations of all B — * W modes simultaneously, 
or are too small to achieve the purpose. When discussing the annihilation effect on penguin- 
dominated decays, we found that the B pK* polarization data can be understood, which 
form the third category mentioned in the Introduction. 

4.1 Annihilation Contribution 

If the power-suppressed annihilation contribution from the {S — P){S + P) penguin operators is 
enhanced by some mechanism, one could have the different counting rules as shown in Eq. (2), 
and it might be possible to reach Rl ~ 0.5. This is the strategy adopted in [14], whose 
analysis was performed in the QCDF approach [15]. In QCDF an annihilation amplitude is 
not calculable due to the end-point singularity, and has to be formulated in terms of several 
free parameters, such as p^. In the current case, different have been introduced for the 
longitudinal and transverse polarization amplitudes in order to fit the data. These parameters 
greatly reduce the predictive power of QCDF. Because varying free parameters to explain 
the data can not be conclusive, we shall estimate the annihilation contribution in the PQCD 
approach, viewing that the PQCD predictions for the penguin annihilation are consistent with 
the measured direct CP asymmetries in B^ K~^7r~, Tr'^n" [11, 12]. Such a calculation for 
a pure-penguin VV mode has been performed in [3]. As shown in Table 4, both the penguin 
annihilation and nonfactorizable contributions help reduce R^. However, the combined effect 
is still not sufficient to lower the fractions Rl of the B (pK* decays down to around 0.5. 

Note that our predicted relative strong phases among Al, A\\, and A± are consistent with 
the B (pK*^ data: 

011 = 2.21 ± 0.22 ± 0.05 (rad.) , (f)± = 2.42 ± 0.21 ± 0.06 (rad) [5] , 

011 = 2.34l°|3 ± 0.05 (rad.) , (f)± = 2.47 ± 0.25 ± 0.05 (rad) [6] . (69) 

Table 4 also implies that R^ can decrease down to 0.75 for a pure-penguin VV mode, after 
taking into account the penguin annihilation and nonfactorizable contributions. Since the 
S+ — > decay is a pure-penguin process, and the p meson mass is not very different from 

the meson mass, the above PQCD analysis applies. Hence, we expect the longitudinal fraction 
Rl ~ 0.75 for the B~^ p^K*^ decay, which is consistent with the Babar measurement, but a 
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bit larger than the Belle measurement. Due to the large uncertainty of the Belle data, there is 
in fact no discrepancy. 



Mode 




I^IlP 




0|| (rad.) 


(f)±{rad.) 




0.923 


0.040 


0.035 


TT 


TT 


(11) 


0.860 


0.072 


0.063 


3.30 


3.33 


(III) 


0.833 


0.089 


0.078 


2.37 


2.34 


(IV) 


0.750 


0.135 


0.115 


2.55 


2.54 


(j)K*+{l) 


0.923 


0.040 


0.035 


TT 


TT 


(11) 


0.860 


0.072 


0.063 


3.30 


3.33 


(III) 


0.830 


0.094 


0.075 


2.37 


2.34 


(IV) 


0.718 


0.133 


0.111 


2.55 


2.51 



Table 4: (I) Without nonfactorizable and annihilation contributions, (II) add only nonfactoriz- 
able contribution, (III) add only annihilation contribution, (IV) add both nonfactorizable and 
annihilation contributions. 

We then come to another mode 5+ — > p^K*^, to which the tree operators contribute. 
Adopt the Babar measurement Rl ~ 0.8 for the — > p^K*^ decay, and assume that the 
tree contribution, which obeys the counting rules in Eq. (1), affects only the longitudinal 
fraction. Hence, we simply add the color- allowed and color-suppressed tree amplitudes T + C ~ 
0.6exp(-90°i)P, which was extracted from the B Ktt data [52], to the B+ p+K*° 
longitudinal polarization amphtude P. The phase —90° has included the weak phase 03 ~ 60°. 
Without an exphcit computation, we derive the polarization fractions for the S+ — > p^K*'^ 
decay, 

Rl ~ 0.86 , R±r^ 0.07 , (70) 

which are consistent with the data within la. We emphasize that we did not attempt a rigorous 
calculation of the B pK* decays here, which deserves a separate paper. In the perturbation 
theories, such as PQCD and QCDF, the factorizable B — > ttK amplitudes and the factorizable 
B — > pK* amplitudes with longitudinal polarizations are very similar. Small differences arise 
only from the meson masses and the distribution amplitudes. Therefore, the estimation using 
the three amplitudes from the B — > ttK modes, which have been available in the literature, 
makes sense. 

The analysis and the result of the modes B — > ujK* should be similar to those of 5 — > p^K* . 
The explicit PQCD analysis of the B — > p{uj)K* polarizations will be performed elsewhere. 
In conclusion, it is not difficult to accommodate the polarization data of the third category, 
the B pK* decays, within the Standard Model by means of the penguin annihilation and 
nonfactorizable contributions. It is also interesting to propose that the measurement of the 
B —>■ LuK* polarizations can test the PQCD approach. 
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4.2 Charming Penguin 

A charming penguin arises from the nonpcrturbative dynamics involved in a charm quark loop 
[53]. It is not calculable, has to be parameterized as a free parameter, and could be as large 
as a leading contribution. Recently, it has been introduced into soft-coUinear effective theory 
(SCET) in order to account for the large B — > 7r°7r° branching ratio [16]. The inputs of the 
measured CP asymmetries St^tt and A^^r demand a complex charming penguin, leading to a 
large penguin-over-trcc ratio \P/T\ ~ 0.7 [16]. With this P/T from the data, a large branching 
ratio B{B^ —>■ 7r°7r°) ~ 1.9 x 10~^ was obtained. SCET does not attempt to explain why \P/T\ 
is so large, even though the Standard Model calculations based on PQCD and QCDF give 
\P/T\ — 0.23-0.29. Another concern is that the S — > tt form factor from the data fitting is as 
small as 0.17 in the presence of the large charming penguin, in conflict with the values 0.28 
from lattice QCD [54] and from light-cone sum rules [55, 56]. 

It has been also proposed that the charming penguin may be large enough to modify the 
counting rules in Eq. (1), and to explain the abnormal B (pK* polarization data [16]. 
However, one also requires different free parameters for the different helicity amplitudes in 
order to lower the longitudinal polarization fraction, and to enhance the transverse polarization 
fractions. In this sense, SCET is similar to QCDF [14], where different pA were introduced for 
the different helicity amplitudes. One needs different parameters for different modes too, such 
as i? — > pK* and B (pK* . Our comment on SCET is then the same as on QCDF: the 
explanation by introducing as many parameters as necessary is always plausible, but can not 
be conclusive. We point out that the current SCET formalism is only of leading power: the 
chirally enhanced terms, proportional to mo/ms, have been dropped, and the annihilation (or 
H^-exchange) amplitudes have not yet been formulated. We speculate that if the annihilation 
amplitude is included into SCET, the charming penguin may not be so essential. On the other 
hand, the charm-loop correction is well-behaved in perturbation theory without any infrared 
singularity, which has been known as the Bander-Silverman-Soni mechanism [57] , implying that 
its nonpcrturbative piece is unlikely to be large. Besides, the light-cone-sum-rule analysis has 
supported a small charming penguin [58]. 

We have taken this chance to investigate the charm-quark loop correction to the B — > (pK* 
polarization fractions in the PQCD approach. The gluon invariant mass attaching the charm- 
quark loop can be defined unambiguously as 

= {l- X2)x^m% - |k2T - ^^t? , (71) 

with X2 and k2T {xz and k^T) being the momentum fraction and the transverse momentum in 
the K* (0) meson, respectively. It turns out that this effect increases i?L by about 5%, and is 
negligible. It also decreases the relative strong phases 4>\\ and 0_l a bit. 

4.3 Rescattering Effect 

It has been proposed to explain the B (j)K* polarization data through the rescattering effect 
[17, 18, 19], 

B D^;^D^*^ (f)K* . (72) 
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The motivation is that the longitudinal polarization fraction of the intermediate states D*D*, 
as low as 0.5, might propagate into the final state (f)K*. First, the massive B meson can decay 
into (f)K* through many intermediate states. The analysis in [17, 18, 19] was restricted to only a 
few channels, and likely to be model-dependent [22]. The truncation of the higher intermediate 
states in this kind of analyses has been criticized [21]. Second, if this mechanism works for the 
B — > (j)K* modes, it will also work for B — > pK*, which involve the same intermediate states. 
As obtained in [19], Rl of both the S+ p'^K*^ and S+ p^K*^ decays are as low as 0.6. 
This observation is expected, since the additional tree amplitudes in the latter can not change 
Rl very much. However, the data in Tabic 1 indicate Rl ~ 0.96 for the B^ p^K*^ decay. In 
other words, the B^ p^K*^ polarization data, obeying the naive coimting rules, have strongly 
constrained the rescattering effect. Third, the D*D and DgD* intermediate states, contributing 
to the P-wave component, could affect the perpendicular polarization of the B — > (j)K* decays. 
Unfortunately, there exists a strong cancellation among these two channels due to the CP and 
SU(3) (CPS) symmetries [19]. The D*D* intermediate state survives the CPS symmetry, which, 
however, exhibits a vanishing R± as in Eq. (17). Therefore, the rescattering effect leads to the 
pattern, 

Rl ~ R\\ > R± , (73) 

contrary to the observed approximate equality R\\ ~ R±. 

Furthermore, it has been known that the B — > KK decays are sensitive to rescattering 
effects. The B — > KK branching ratios measured recently well agree with the PQCD predictions 
[59] as shown in Table 5, leaving very limited room for the rescattering effect. Note that no 
theoretical errors were presented in [59], since the detailed investigation of uncertainties in the 
PQCD approach was available only after Ref. [45]. Roughly speaking, the theoretical errors on 
PQCD predictions for branching ratios of two-body charmless B meson decays are about 30%. 
Viewing the contradiction of Eq. (73) to the B — > (f)K* polarization data, and the constraints 
from the measured B~^ — > p^K*^ polarizations and from the mcasiucd B — >■ KK branching 
ratios, we intend to conclude that the rescattering effect is not a satisfactory resolution to the 
polarization puzzle. 



Branching Ratio 


PQCD 




Babar [60] 


B{B+ K+K"") 
B^B'' K^K'') 


1.65 X IQ-^ 
1.75 X 10-6 


(1.45: 
(1.19: 


^''oii±0.11) X 10-6 
±0.13) X 10-6 



Table 5: PQCD predictions for the CP-averaged B — > KK branching ratios and the data. 



4.4 Magnetic Penguin 

If the B — > pK* data can be understood in the Standard Model by means of the penguin 
annihilation and nonfactorizable contributions, and only the B — > (f)K* decays exhibit an 
anomaly, it is natural to look for a unique mechanism for the latter. Such a mechanism, the 
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Figure 4: Some diagrams from the magnetic penguin which contribute to the transverse polar- 
ization fraction of the B (pK* decays. 



b ^ sg transition, has been proposed in [20]. The novel idea is that the transversely polarized 
gluon from the transition propagates into the meson, enhancing the transverse polarization 
amplitudes. The relevant matrix element was then parameterized in terms of a dimensionless 
free parameter k. Assuming this parameter to be k ~ —0.25, the authors of [20] claimed 
that the B (pK* polarization data could be accommodated within the Standard Model. 
Similarly, varying a free parameter to fit the data can not be conclusive, and a reliable estimate 
of the K value is necessary. As pointed out in [20], the same mechanism also contributes to the 
B —>■ ujK* decays, and small Rl ~ 0.5 have been predicted. Therefore, the measurement of the 
B ujK* polarizations will impose a stringent test on this proposal in the future. Note that 
PQCD postulates, contrary to [20], that Rl of the B ujK* decays are as large as those of 
B p^K*. 

Besides the above experimental discrimination, we shall estimate the order of magnitude 
of K in the framework of FA, following the method in [61]. The weak effective Hamiltonian 
contains the b ^ sg transition, 

-^V:,V,,Cs,Os, , (74) 

with the magnetic penguin operator, 

Osg = £^mks,a^,{l + j,)^G''^%, , (75) 

i, j being the color indices. The picture described in [20] is displayed in Fig. 4: one or more 
coUinear gluons, emitted from the B K* form factor, produce the s and s quarks in the color- 
octet state. They, together with the transversely polarized gluon from the b ^ sg transition, 
fragment into the color-singlet transversely polarized meson. According to this picture, 
we introduce three-parton distribution amplitudes to absorb the nonperturbative dynamics 
associated with the (p meson. Another diagram, in which the transversely polarized gluon 
produces the s and s quarks and the collinear gluon flows into the meson directly, does not 
contribute, since the transverse polarization of the meson is mainly carried by its gluonic 
parton. 

We first argue that the leading picture in Fig. 4(a), where the collinear gluon attaches 
the s quark from the b —>■ sg transition, diminishes due to the G-parity. The corresponding 



20 



three-parton distribution amplitudes are defined via the matrix elements [32], 



{<l>{Ps,eliT))\s{-z)gG,,{vzMz)\0) 



= -iP: 



3a 



P^AuiT) - el,{T)Ps, flV{v, Ps-z), 



{<P{Ps,el{T))\s{-z)gG^,{vz)j^j,siz)\0) 



P 



3a 



Ps^JT) - Ps,el{T) /3^.i(t;, P, ■ z) , 



(76) 



(77) 



with the dual gluon field strength tensor (5^,^ = e^^pfjC"^ /2. Other three-parton distribution 
amplitudes, irrelevant to the discussion below, are not quoted here. The Fourier transformation 
of the distribution amplitude V gives 



n.P3..)^/He.p[,.3..(..-.....,M....,..,). 



(78) 



with Xg, Xg, and Xg being the momentum fractions carried by the s quark, the s quark, and the 
gluon, respectively, and the integration measure. 



[[dx] = [ dxg I dxg [ 
J Jo Jo Jo 



s 1^ dxg6 ( 1 



E 



Xi 



(79) 



The Fourier transformation of A is defined in the same way. The three-parton twist-3 vector 
meson distribution amplitudes have been also studied in [62] . The asymptotic models of V and 
A have been parameterized as [32, 63, 64] 



V{xs, Xs, Xg) = 5040(x^ - 



'^S J'^ S'^ S'^ g 1 



(80) 
(81) 



with the shape parameter ujf-Q = —2.1. The antisymmetry (symmetry) of V (A) between the 
exchange of Xs and Xs is a consequence of the G-parity transformation in the SU(3) limit [32]. 

fV 
'30 



The constant f^^ is chosen such that V is normahzed according to 



y [dx] (Xs - Xs)V(Xs, Xs, Xg) 



(82) 



We factorize the matrix element in Eq. (76) and the B K* transition form factor out of 
Fig. 4(a). The remaining part is the hard kernel, which must be symmetric under the exchange 
of Xs and Xs- Therefore, Fig. 4(a), written as the convolution of the symmetric hard kernel 
with the antisymmetric three-parton distribution amplitude V , vanishes. The diagram with 
the coUinear gluon attaching the h quark does not contribute, because its hard kernel is also 
symmetric. No matter how many infrared gluons are involved in the s-s quark pair production, 
there is no contribution for the same reason. If factorizing the matrix element in Eq. (77) out of 
Fig. 4(a), it is easy to find that the corresponding hard kernel vanishes, because the s-s quark 
pair does not form an axial-vector current. Therefore, we conclude that the leading diagram 
does not contribute to the transverse polarization amplitudes of the B (f)K* decays. 
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To survive the above suppressions, one has to consider subleading diagrams such as Fig. 4(b), 
in which both the s quark and the gluon from b ^ sg flow into the (j) meson, or such as Fig. 4(c), 
in which one more infrared gluon fragments into the (f) meson. For Fig. 4(b), an extra hard gluon 
is necessary for producing the s-s quark pair, such that the price to pay is the suppression. 
For Fig. 4(c), four-parton distribution amplitudes are involved, whose contribution is power- 
suppressed. We shall show that the order of magnitude of k from Fig. 4(b) is, unfortunately, 
as small as 0.01, far away from k ~ —0.25 required by the data. For simplicity, we analyze the 
contribution from the three-parton distribution amplitude V, and our conclusion applies to that 
from A. Due to the lack of the information of the four-parton twist-4 distribution amplitudes, 
we can not estimate the contribution from Fig. 4(c) in a reliable way. However, it is of higher 
power in m^/ms, and unlikely to be huge. 

Insert the Fierz identity. 



lijllk = ^(7a)ife(7")/j H , 

where the irrelevant terms have been suppressed, and the identity for color matrices. 



(83) 



(84) 



to change the fermion and color flows of the outgoing s and s quarks, respectively. Figure 
4(b), where the additional hard gluon attaches the s quark going into the K* meson, is then 
factorized into 



Gi 
1 

x-tr 
4 



M = --^VZVt,^m,CsgJ[dx]lFT{<p{P,,e;{T)M^^^^^ 

2tr{T^T^T^T^) 



•■■77 ixjT, — rV\i^M^(i + 75) 



(P2 + xM^ ' 



(85) 



where IFT means the inverse Fourier transformation. The indices A denote the hard gluon 

vertices, and the dots in the trace represent the Feynman rules associated with the B — > K* form 
factor. Employing 2tr{T^T^T^T'') = -5''7(2A^c), Eq. (76) and Eq. (78), the above expression 
becomes 
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x-tr 
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M = -^VtlVtb:^^^mbCsgflf, J[dx]V{xs,Xs,Xg 

;^crnv{l + 75) ' 



•7 /P37A 
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(86) 



Neglecting the light meson masses, assuming nib 
Dirac matrices in the trace, we derive 



niB, and working out the product of the 



M = 



G 



a. 



x{K*-(P,,el(T))\sia,4l + ^,)b\B-(P,))P^er(T). 



(87) 
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Comparing the above expression with Eq. (4) in [20], the parameter k is given by 

with = /|^/ {f^m,^. For the values = 0.4 and C]^ = 0.013 [32], and the model distribution 
amplitude in Eq. (80), wc obtain k ^ 0.004. Other diagrams with the hard gluon attaching the 
b quark and the transversely polarized gluon can be analyzed in a similar way, and the results 
are of the same order of magnitude. Adding these contributions leads to 

K ^ 0.01 . (89) 

Hence, we intend to conclude that the magnetic penguin is not sufficient to resolve the B — > (f)K* 
puzzle. 



5 CONCLUSION 

In this paper we have investigated most of the B VV modes carefully. Our observation is 
that the B — > VV modes can be classified into four categories. For the tree-dominated decays, 
the polarization fractions are basically determined by kinematics. The upper part of Table 2 can 
be understood by kinematics in the heavy-quark limit. The longitudinal polarization fractions 
follow the mass hierarchy among the D*, D* and p (K*) mesons. The lower part of Table 2 can 
be understood by kinematics in the large-energy limit. We always have Rl ~ 1 for the decays 
into two light vector mesons. It has been found that the above simple kinematic estimates 
in the heavy-quark and large-energy limits are robust under subleading corrections. For this 
part, we have analyzed the next-to-leading-power corrections to the universal IW function, 
the nonfactorizable contributions, the two-parton twist-4 contributions, and part of next-to- 
leading-order contributions from the charm-quark loop, all of which are negligible. That is, 
QCD dynamics plays only a minor role for the polarizations of the tree-dominated decays. 

As a byproduct, we have predicted the longitudinal polarization fractions of the S+ — > 
(D*"*", D*"'")p° modes in the large-energy limit using FA, and found Rl ~ 0.7. We have also 
calculated the polarization fractions of the B'^ jj*Ojj*o ^jg^ay explicitly in the PQCD ap- 
proach based on factorization theorem. The result Rl ~ 1 is quite different from Rl of 
other B — > D*D* modes, since it is dominated by the nonfactorizable ly-exchange topology. 
The above predictions can be confronted with the future data. 

For the penguin-dominated modes, the polarization fractions can deviate from the naive 
counting rules based on kinematics, because of the important annihilation contribution from 
the {S — P){S + P) operators. This mechanism explains the third category listed in the 
lower part of Table 1: Rl can decrease to 0.75 for the B^ — > p'^K*^ mode. Adding the tree 
contribution, Rl of the S+ — > p^K*^ decay can go up to about 0.9. We have postulated from 
the viewpoint of PQCD that the B — > ujK* decays also belong to the third category, and should 
show Rl similar to those oi B ^ p^K* . All the above three categories can be accommodated 
within the Standard Model. Only the fourth category, the B (pK* decays, can not. They are 
dominated by the penguin contribution, but their Rl ~ 0.5 are much lower than 0.75. We have 
carefully analyzed the various mechanism proposed in the literature to resolve this anomaly. 
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and concluded that none of them is satisfactory. Therefore, the B (pK* polarization data 
remain as a puzzle. However, we emphasize that we are not claiming a signal of new physics, 
since the complicated QCD dynamics in the B — > VV decays has not yet been fully explored. 
For example, a smaller B ^ K* form factor Aq could decrease Rl significantly [65]. 
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